function res = cal_CDcoes2(p)
res=[1];
if p == 4
    res = [res, -1.0/12];
end
    